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Introduction of a Josephson field effect transistor (JoFET) concept [I] sparked ac¬ 
tive research on proximity effects in semiconductors. Induced superconductivity and 
electrostatic control of critical current has been demonstrated in two-dimensional 
gases in In As pi [3], graphene [4] and topological insulators |5|? -[8], and in one¬ 

dimensional systems [9 - FTl ] including quantum spin Hall edges [12, 13]. Recently, 
interest in superconductor-semiconductor interfaces was renewed by the search for 
Majorana fermions [14 > 815] , which were predicted to reside at the interface[16 
More exotic non-Abelian excitations, such as parafermions (fractional Majorana 
fermions) [TI)-21] or Fibonacci fermions may be formed when fractional quantum 
Hall edge states interface with superconductivity. In this paper we develop transpar¬ 
ent superconducting contacts to high mobility two-dimensional electron gas (2DEG) 
in GaAs and demonstrate induced superconductivity across several microns. Super¬ 
current in a ballistic junction has been observed across 0.6 ^m of 2DEG, a regime 
previously achieved only in point contacts but essential to the formation of well sep¬ 
arated non-Abelian states. High critical fields (> 16 Tesla) in NbN contacts enables 
investigation of a long-sought regime of an interplay between superconductivity and 
strongly correlated states in a 2DEG at high magnetic fields p2K2Ti |. 

Proximity effects in GaAs quantum wells have been intensively investigated in the past and 
Andreev reflection has been observed by several groups [28H31] . Unlike in InAs, where Fermi level 
(Ep) at the surface resides in the conduction band, in GaAs Ep is pinned in the middle of the 
gap which results in a high Schottky barrier between a 2DEG and a superconductor and low 
transparency non-ohmic contacts. Heavy doping can move Ep into the conduction band and, 
indeed, superconductivity has been induced in heavily-doped bulk n ++ GaAs[32j. In quantum 
wells similar results were obtained by annealing indium contacts [33], however the critical field of 
indium is ^ 30 mT which is well below the fields where quantum Hall effect is observed. 

In conventional quantum well structures AlGaAs barrier between 2D electron gas (2DEG) and 
the surface of the sample adds an extra 0.3 eV to the Schottky barrier when contacts are defused 
from the top. We alleviated these problems by growing an inverted heterojunction structures, 
where a 2DEG resides at the GaAs/AlGaAs interface but the AlGaAs barrier with modulation 
doping is placed below the 2DEG, see Fig. [lj Contacts are recessed into the top GaAs layer in 
order to bring superconductor closer to the 2DEG. A thin layer of AuGe and NbN superconductor 
form low resistance ohmic contacts to the 2DEG after annealing. The inverted heterostructure 
increases the contact area of side contacts compared to quantum well structures by utilizing all 
GaAs layer above the heterointerface for carrier injection (130 nm in our inverted heterostructure 
vs 20 — 30 nm in typical quantum wells). 

Here we report induced superconductivity in two devices from different wafers, sample A has 
long (70 /im) contacts separated by 1.6 jim of 2DEG, for sample B contacts are formed to the edge 
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FIG. 1. Devices design and superconducting transition, (a) Scanning electron microscope (SEM) 
images of test devices similar to samples A and B. Enlarged region for sample B is an atomic force microscope 
(AFM) image of a real sample. 2D gas regions are false-color coded with green, superconducting and normal 
contacts are coded with orange and blue, respectively, (b) Simulation of the conduction band energy profile 
in the heterostructure [34l? ]. (c) T-dependence of resistance between contact 3 and 4 in Sample B measured 
with 10 nA ac excitation. Superconducting transition is observed at T c « 290 mK. 


of a mesa with 0.6 gm separation. Details of device fabrication are described in Methods. When 
cooled down to 4 K in the dark both samples show resistance in excess of 1 MD. After illumination 
with red light emitting diode (LED) a 2DEG is formed and 2-terminal resistance drops to < 500D. 
As shown in Fig. |TJi sample resistance gradually decreases upon cooldown from 4 K to the 

base temperature and the S-2DEG-S junctions becomes superconducting at T c ~ 0.3 K. 

Voltage-current V(I) characteristics for two S-2DEG-S junctions (between 8-9 for sample A and 
between 3-4 for sample B) are shown in Fig. [2j Both samples show zero resistance state at small 
currents with abrupt switching into resistive state at critical currents I c — 0.22 /jl A and 0.23 fi A 
for samples A and B respectively. V(I) characteristics are hysteretic most likely due to the Joule 
heating in the normal state. 

The most attractive property of a high mobility 2DEG is large mean free path l £o, with l = 
24 fim and the BCS coherence length £o = hvf/irA = 0.72 gm for sample B. Here Vf = hy/27rn/m is 
the Fermi velocity, n is a 2D gas density, m is an effective mass, and A = 1 .T^ksTc = 46 /ieV is the 
induced superconducting gap. Evolution of V (/) with T is shown Fig. [3^. Experimentally obtained 
T-dependence of I c is best described by the Kulik-Omelyanchuk theory for ballistic junctions 
(L <C l)\ 35], the blue curve on Fig. If- For comparison we also plot 7 C (T) dependence for the dirty 
limit L <C VWo [36] . which exhibits characteristic saturation of I c at low temperatures. 

In short ballistic junctions L <C £o l the product I c (0)Rn = 7rA/e does not depend on the 
junction length L. For L ~ £ 0 this product is reduced by a factor 2£o/(T + 2£o) [ 32 ] • The measured 
I c Rn — 83 /iV for sample B is in a good agreement with an estimate 7rA/e*2^o/(T + 2^o) = 90 /iV. 
For sample A the I c Rn = 19 /iV while the estimated product is ~ 50 /iV. The reduction is consistent 
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FIG. 2. Induced superconductivity in a high mobility 2D electron gas in GaAs. Voltage-current 
characteristics and differential resistance are measured between 8-9 for sample A and between 3-4 for sample 
B at base temperature , dV/dl is measured with I ac = 1 nA. Induced superconductivity with zero voltage 
is observed with critical currents I c ~ 220 nA for sample A and I c ~ 230 for sample B. 


with the geometry of sample A, where a region of the 2DEG with induced superconductivity is 
shunted by a large region of a 2DEG in a normal state. 

Transparency of superconducting contacts can be estimated from the suppression of the su¬ 
perconducting gap in the S-2DEG-S junction between 3-4 in sample B. In one-dimensional junc¬ 
tions the induced gap A = Aq r ^ Ao depends on the broadening of Andreev levels within the 

semiconductor [ 3jjj T = where we introduce contacts transparencies D\ and D 2 . We 

assume for simplicity that D\ — D 2 = 1/(1 + Z 2 \ where 0 < Z < 00 is a interface barrier strength 
introduced in [38], and Bagwell’s effective channel length L e ff = L + 2£o- Using NbN supercon¬ 
ducting gap Aq = 2.02fcgTjS, (NbN is a strong-coupling superconductor, T ^ = 11 K) we obtain 
Z = 0.2. This value is consistent with the fit of the I c vs T dependence with D as a free parameter, 
see Supplementary Material for details. Similar values of Z can be estimated from the analysis of 
the shape of dI/dV(V) characteristics at elevated temperatures, as shown in Fig. [3} At T < T® 
Andreev reflection at S-2DEG interfaces results in an excess current flowing through the junction 
for voltage biases within the superconducting gap Ao/e and corresponding reduction of a differ¬ 
ential resistance dV/dl by a factor of 2. In the presence of a tunneling barrier normal reflection 
competes with Andreev reflection and reduced excess current near zero bias, resulting in a peak 
in differential resistance. Within the BTK theory[38] a flat dV/dI(V) within Ao/e, observed in 
our experiments, is expected only for contacts with very high transparency Z < 0.2. For larger 
Z > 0.2 a peak at low biases is expected (see Supplementary Material). Several features of the 
experimental I(V) need to be mentioned. First, we observe several sharp peaks in the resistance 
at high biases (around 2 mV and 4 mV for T = 4K). Similar sharp resonances has been observed 
previously @ni, where authors attributed their appearance to the formation of Fabry-Perot reso¬ 
nances between superconducting contacts. In our devices the superconducting region is shunted by 
a low resistance (< 10012) 2DEG, thus appearance of > 10 k 12 resonances cannot be explained by 
resonant electron trapping between contacts. These resonances are also observed in I(V) character¬ 
istics of a single S-2DEG interface (measured in the S-2DEG-N configuration between contacts 3-6, 
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FIG. 3. Temperature dependence of superconductivity in a ballistic junction, (a) Evolution of the 
induced superconductivity with T for sample B. The R(I) curves are offset proportional to T for T > 50 mK. 
(b) Temperature dependence of critical current J C (T) is extracted from (a) and compared to the expected 
T-dependence for ballistic and diffusive regimes (reduced I c compared to Fig. [2] is due to larger I ac = 10 nA 
used in this experiment), (c) High-temperature data shows Andreev reflection (excess current and reduced 
dV/dl around V = 0. The curves are not offset. In (d) excess current is modeled within the BTK theory[38] 
with Z = 0.2. 


see Supplementary Material Fig. S2). Differential resistance does not change substantially across 
resonances, ruling out transport through a localized state. We speculate that in the contacts where 
these resonances are observed superconductivity is carried out by quasi-ID channels, and jumps 
in I/V characteristics are due to flux trapping at high currents. This scenario is consistent with 
the observation that peaks shift to lower currents at higher fields, see Fig. [4j The second notable 
feature of our data is reduction of the zero-bias resistance by ~ 2.6 at low temperatures, while 
Andreev reflection limits the reduction to the factor of 2. We attribute this reduction to the mul¬ 
tiple Andreev reflection between two closely-spaced contacts, for contacts with larger separation 
( 20 fi m) multiple Andreev reflection is suppressed and the reduction of resistance by a factor of 2 
is observed, see Supplementary Material Fig. [S2j 


Finally, we present magnetic field dependence of induced superconductivity. The low-field data 
is shown in Fig. |4](a,b), where black regions correspond to zero differential resistance. Induced 
superconductivity is suppressed at « 0.2 T in both samples. In sample A a narrow region of a 
2DEG with induced superconductivity is confined between large NbN superconducting leads with 
rigid phases. Perpendicular magnetic field twists the phase in the 2DEG resulting in Fraunhofer- 
like oscillations of the critical current. In this sample, though, the 2DEG extends beyond the 
narrow region between the contacts and I c does not decrease to zero and abrupt jumps in I c 
reflect multiple flux jumps. The period of oscillations is ~ 0.5 mT which corresponds to an area 
of 4.1 gm 2 , much smaller than the area of the 2DEG between the contacts (« 120 gm 2 ). This 
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FIG. 4. Magnetic field dependence of induced superconductivity. (a,b) Differential resistance is 
measured as a function of B and Id c for two samples at 40 mK. Induced superconductivity (black region) 
is observed up to 0.2 Tesla in both sample, (c) 3-terminal resistance for a sample with all normal contacts 
(red) and between normal and superconducting contacts in sample B [I (2 — 4) and V (4 — 1) in Fig. [I] is 
measured at 70 mK and 40 mK respectively. B < 0 (B > 0) induces clockwise (counterclockwise) chiral 
edge channels, note resistance scales difference for two field directions. 


observation is consistent with the reduced I c Rn product measured for this sample as discussed 
above. In sample B contacts are fabricated along the edge of the mesa and 2D gas is not enclosed 
between the contacts. Consequently, I c is a smooth function of B. 

Competition between superconductivity and chiral quantum Hall edge states is shown in Fig. [4]^, 
where resistance is measured in a 3-terminal configuration over a wide range of magnetic fields. 
Simple Landauer-Buttiker model of edge states predicts zero resistance for negative and quantized 
Hall resistance for positive field direction for IQHE and FQHE states, which is clearly seen in a 
sample with all normal ohmic contacts (red curve). When a superconducting contact serves as a 
current injector (blue curve), integer v — 1 and fractional v — 2/3 and 3/5 states are well developed 
for B < 0, while the same states are not quantized at proper QHE values for B > 0. If we assume 
that current injection via a superconducting contact results in an extra voltage offset at the contact 
V Q ff « A i n d/e, the measured voltage will be reduces by V Q ff. The magenta bars for B > 0 indicate 
corrected resistance {V — V Q ff)/I for V Q ff = 140 /iV. While this offset may explain the measured 
values for fractional states, a twice smaller V Q ff is needed to reconcile the resistance at v — 1. 
Note that induced gap is smaller at higher B. At low fields states v — 3, 4 and 5 have resistance 
minima for B < 0 indicating partial equilibration of chiral edge currents with the superconducting 
contact, while resistance near v — 2 has a maximum. Zero resistance at v — 1 and large resistance 
at v — 2 are in contrast to the theoretical prediction that v — 2 state should be stronger coupled 
to a superconducting contact than v — 1 [23] . 

Methods 

The GaAs/AlGaAs inverted heterojunctions were grown by molecular beam epitaxy (MBE) on 
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semi-insulating (100) GaAs substrates with the heterointerface placed 130 nm below the surface 
and 5-doping layer 30-40 nm below the GaAs/AlGaAs interface. Samples were fabricated from 
two wafers with density and mobility n — 2.7 x 10 11 cm -2 , fi — 2 x 10 6 V-s/cm 2 (sample A) 
and n = 1.7 x 10 11 cm -2 , /i = 4 x 10 6 V-s/cm 2 (sample B). Superconducting contacts were 
defined by standard electron beam lithography. First, a 120 nm - deep trench was created by 
wet etching. Next, samples were dipped into HCkH^O (1 : 6) solution for 2 s and loaded into a 
thermal evaporation chamber, where Ti/AuGe (5nm/50nm) was deposited. Finally, 70 nm of NbN 
was deposited by DC magnetron sputtering in Ar/N 2 (85%/15%) plasma at a total pressure of 2 
mTorr. The deposition conditions were optimized for producing high quality NbN films (T c = 11 
K and B c > 15 Tesla) with minimal strain[4T]. The contacts were annealed at 500° C for 10 min in 
a forming gas. The measurements were performed in a dilution refrigerator with base temperature 
< 30 mK, high temperature data was obtained in a variable temperature 3 He system. Samples 
were illuminated with red LED at 4 K in order to form a 2D gas, 2-terminal resistance drops from 
> 1MD before illumination to < 500 D after illumination. 
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TEMPERATURE DEPENDENCE OF THE CRITICAL CURRENT 
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FIG. SI. Analysis of the temperature dependence of the critical current. Scaled (a) and unsealed 
(b) product I c Rn is calculated using Eq. (SI) for different transparencies D and a = 1. Red dots are 


experimental data. Dashed line in (b) is for a = 0.7 and D = 1. In (c) root-mean-square deviation between 
the best fit and the experimental data is shown for different D, coherence length £ obtained from the best 
fit are red triangles. 


Haberkorn et al. m generalized Kulik-Omelyanchuk current-phase relations [35 \ 36j to the case 
of arbitrary transparency of a tunnel barrier D inserted into the Josephson junction by directly 
solving Gor’kov’s equations. They obtain the following current-phase relation: 
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where A(T) is the BCS gap. For a = 1 this equation interpolated between diffusive (D — 0) and 
ballistic (D — 1) junctions. Critical current can be found as I c {T)Rjsf = max[I s {(j)^ T)Rn], We 
introduce coefficient a to account for the reduction of the critical current due to the finite length of 
the junction L, a — 2£/(L+2£) [57]. The best fit of the experimental I c Rjsr(T) dependence assuming 
both a and D as free parameters is obtained for D = 1 and a = 0.7, see Fig.[Sl|a,b). For the contact 
spacing L — 0.63 /m i this a corresponds to £ = 0.76 /im, consistent with the BCS coherence length 
£o = hvf/nA = 0.72 gm. Transparency D can be related to the dimensionless barrier strength Z 
introduced in the Blonder-Tinkham-Klapwijk (BTK) theory [38]. D — 1/(1 + Z 2 ), and the fit sets 
the upper limit on Z, Z < 0.1. The quality of the fit parameters can be assessed from Fig. [SlYc) , 
where rms error for the best fit with a fixed D and a as a free parameter (rms deviation) 2 = 
J2i{[I c (Ti)RN] theory — [I c (Ti)R]y ] exp } 2 is plotted for different D. The rms deviation has a clear 
global minimum at D 1. Note that the coherence length for D < i, obtained from the fitting 
parameter cq becomes smaller than the estimated £q. 
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FIG. S2. Temperature dependence of differential resistance. Left 6 plots: normalized differential 
resistance is calculated using BKT theory, Eq. [S2] for different barriers Z and temperatures between 4 
and 11 K with a step of 1 K. Right 2 plots: experimentally measured differential resistance between two 
superconducting contacts (R 3 - 5 ) and a normal-superconducting contact (R 4 - 7 ) in sample B (the normal 
contact has high resistance). 


ANALYSIS OF EXCESS CURRENT ABOVE THE INDUCED SUPERCONDUCTIVITY 

GAP 


Transparency of the superconductor-semiconductor interface can be estimated from the shape 
of the dV/dI(V ) characteristic, where competition between Andreev and normal reflections results 
in a peak in differential resistance when a tunneling barrier is present at the superconductor- 
semiconductor interface (transmission D = 1/(1 + Z 2 ) < 1). Differential resistance for different 
temperatures can be calculated using Blonder-Tinkham-Klapwijk (BTK) theory[38]: 


dl_ 

dV 



df 0 (E-eV) 

d{eV) 


[1 + A(E)-B(E)]dE, 


(S2) 


where fo(E) is the Fermi Dirac function and A{E) and B(E) are energy-dependent Andreev and 
normal reflection coefficients, respectively. Both coefficients depend on the gap of NbN Ao = A(T) 
with — 11 K and the interface barrier strength Z. In Fig. S2 we plot differential resistance 
for different values of Z. At low T for Z = 0 the barrier is transparent (D = 1) and all incident 
electrons are Andreev reflected, which leads to the a reduction of differential resistance by a factor 
of 2 within the energy gap Aq. When Z is finite, part of the incident electrons undergoes normal 
reflection which results in the increase of the resistance within the gap. 

The exact shape of experimental curves differ from the shape predicted by the BKT theory, the 
most important deviation being sharp minima near V = 0 observed at T close to Tj? as compared 
to a much smoother BKT dependence. To account for a similar sharpening of a zero-bias peak 
in less transparent contacts (Z > 2) it has been assumed that a thin normal region is formed 
between NbN contacts and a 2DEG[43J. This more elaborate theory introduces two more fitting 
parameters for the superconducting-normal and normal-2DEG interfaces, but does not change the 
main qualitative prediction of a simpler BTK theory: appearance of a peak near V = 0 for Z > 0.2 
in dV/dI(V ) characteristics. 
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Experimentally, we observe no zero-bias peak in dV/dl(V) characteristics measured between 
two superconducting contacts R 3-4 (S-2DEG-S) or between superconducting and normal contacts 
i? 8_9 (S-2DEG-N), see Fig. [3] and [S2| thus we can set an upper limit Z < 0.2 and lower limit 
D > 0.96 for our contacts. 


COMPARISON BETWEEN CONVENTIONAL HETEROSTRUCTURES AND 
INVERTED SINGLE INTERFACE HETEROSTRUCTURE 


Comparison between conventional heterostructures and inverted heterostructures used in this 
work is shown in Fig. S3 In conventional quantum well (b,d) and single interface heterostructures 
AlGaAs barrier between 2D gas and the surface adds 0.3 eV to the Schottky barrier if contacts 
are defused from the surface. For side contacts inverted single heterointerface (a) increases the 
exposed GaAs cross section for Cooper pair injection. 
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FIG. S3. Comparison between conventional heterostructure and inverted single interface het¬ 
erojunction. Conduction band profile is plotted for (a) inverted single interface heterojunction used in our 
experiments and typical (b) modulation-doped quantum well, (c) single heteroj unction, and (d) inverted 
quantum well. Dash lines indicate position of modulation doping. 



































